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Using the density functional theory (DFT) with the relativistic mean field (RMF) 
' model, we study the non-uniform state of nuclear matter, "nuclear pasta". We 

self-consistently include the Coulomb interaction together with other interactions. 
It is found that the Coulomb screening effect is significant for each pasta structure 
but not for the bulk equation of state (EOS) of the nuclear pasta phase. 



1. Introduction 

One of the most interesting features of low-density nuclear matter is the pos- 
sibility of the existence of non- uniform structures, called "nuclear pastas".^ 
At low densities, nuclei in matter are expected to form the Coulomb lat- 
tice embedded in the neutron-electron seas, so as to minimize the Coulomb 
interaction energy. On the other hand, another possibility has been dis- 
cussed: the stable nuclear shape may change from sphere to cylinder, slab, 
cylindrical hole, and to spherical hole with increase of the matter density, 
and "pastas" are eventually dissolved into uniform matter at a certain nu- 
cleon density close to the saturation density, ps ~ 0.16 fm~^. The existence 
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of such "pasta" phases, instead of the ordinary crystaUine lattice of nuclei, 
would modify several important processes in supernova explosions by chang- 
ing the hydrodynamic properties and the neutrino opacity in the supernova 
matter. Also expected is the influence of the "pasta" phases on star quakes 
of neutron stars and pulsar glitches via the change of mechanical properties 
of the crust matter. 

Several authors have investigated the low-density nuclear matter using 
various models. Roughly speaking, the favorable nuclear shape is deter- 
mined by a balance between the surface and the Coulomb energies, as has 
been shown by previous studies, where the rearrangement effect on the den- 
sity profile of the charged particles, especially electrons, by the Coulomb 
interaction is discarded. However, the proper treatment of the Coulomb 
interaction should be very important, as it is demonstrated in Ref. 9; the 
screening of the Coulomb interaction by the charged particles may give a 
large effect on the stability of the geometrical structures. 

We have been recently exploring the effect of the Coulomb screening in 
the context of the structured mixed phases in various first order phase tran- 
sitions such as hadron-quark deconfinement transition, kaon condensation 
and liquid-gas transition in nuclear matter. We treat the nuclear "pasta" 
phases as a part of our project, since they can be considered as structured 
mixed phases during the liquid-gas transition in nuclear matter. 

Our aim here is to study the nuclear "pasta" structures by means of a 
mean field model, which includes the Coulomb interaction in a proper way, 
and we figure out how the Coulomb screening effect modifies the previous 
results without it. 

2. Density Functional Theory with the Relativistic 
Mean-field Model 

To study the non-uniform nuclear matter, we follow the density functional 
theory (DFT) with the relativistic mean field (RMF) model, The RMF 
model with fields of mesons and baryons is rather simple for numerical 
calculations, but realistic enough to reproduce main nuclear matter prop- 
erties. In our framework, the Coulomb interaction is properly included in 
the equations of motion for nucleons, electrons and the meson mean fields, 
and we solve the Poisson equation for the Coulomb potential Vcoui self- 
consistently with them. Thus the baryon and electron density profiles, as 
well as the meson mean fields, are determined in a way fully consistent with 
the Coulomb potential. Note that our framework can be easily extended to 
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other situations. For example, if we take into account meson condensations, 
which are hkely reaUzed in a high-density region, we should only add the 
relevant meson field terms. In Ref. 11 we have included the kaon degree of 
freedom in such a treatment to discuss kaon condensation in high density 
regime. 

To begin with, we present the thermodynamic potential of the form. 
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where = - Me + Vboul - QwN'^O - gpNPO, Vn = /J-B - QloNOJo + 

QpNPo, m*B = TUB - gcrN<J, U{a) = \hmB{gaN<y)^ + \c{gc,NcrY- Here 
we used the local-density approximation for nuclcons and electrons, while 
one still should carefully check its validity. The introduction of the density 
variable is meaningful, if the typical length of the nucleon density varia- 
tion inside the structure is larger than the intcr-nuclcon distance, which 
we assume to be fulfilled. We must also keep in mind that the approxima- 
tion is broken down for small structure sizes, since quantum eff'ects become 
prominent there. For the sake of simplicity we also omitted nucleon and 
electron density derivative terms. If the nucleon length scale were shorter 
than lengths of changes of the meson mean fields, one could not introduce 
the derivatives of the nucleon density but could simplify the problem intro- 
ducing the corresponding contribution to the surface tension. In the given 
case (when we suppressed the derivative terms mentioned above) the result- 
ing nucleon density follows the changes of the meson mean fields. However, 
even in this case the presence of the derivative terms (of the same order 
as for other fields) could affect the numerical results. Here we consider 
large-size pasta structures and simply discard the density variation effect, 
as a first-step calculation, while it can be incorporated in the quasi classi- 
cal way by the derivative expansion within the density functional theory. 
The parameters are set to reproduce the nuclear-matter saturation prop- 
erties. From the variational principle, j^rp: = ^ {<Pi — PO; '^O) ^Coui) or 
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= 0{i = n,p, e), we get the coupled equations of motion as 

- VV + mla = -^+ gMpi''^ + 4'^) (5) 

-V^wo + m^wo = 9u>n{Pp + Pn) (6) 

-V'^Pn + m^pPa= 9pN{Pp- Pn) (7) 
V^Vcoui = 47re^Pch (charge density pch = Pp + Pe) (8) 

P'n = P'B = \J fc|„ + m*g^ + g^Ni^o - QpNPO (9) 

IJ-p = P-B - P-e = \/k^Fp + ™B^ + ffi^iV^O + QpNPO - Vcoul (10) 

Pe = -(Me - V^CouOVStt'- (H) 



Note that first, the Poisson equation (8) is a highly nonhnear equation for 
Vcoui, since pch in r.h.s. includes it in a complicated way, and secondly, the 
Coulomb potential always enters the equation through the gauge invariant 
combination, /ig — Vcoui- 

To solve the above coupled equations numerically, we use the Wigner- 
Seitz cell approximation: the space is divided into equivalent cells with 
some geometrical symmetry. The shape of the cell changes: sphere in 
three dimensional (3D) calculation, cylinder in 2D and slab in ID, respec- 
tively. Each cell is charge-neutral and all the physical quantities in a cell 
are smoothly connected to those of the neighbor cell with zero gradients at 
the boundary. The cell is divided into grid points (iVgrid ~ 100) and the 
differential equations for fields are solved by the relaxation method with 
constraints of given baryon number and charge neutrality. 

3. Bulk Property of Finite Nuclei 

Before applying our framework to the problem of the pasta phases in nuclear 
matter, we check how it can describe finite nuclei. In this calculation, the 
electron density is set to be zero and the boundary condition or the charge- 
neutrality condition is not imposed. However, we assumed the spherical 
shape of nuclei. In Fig. 1 (left panel) we show the density profiles of some 
typical nuclei. We can see how well our framework may describe the density 
profiles of these nuclei. To get a better fit, especially around the surface re- 
gion, we might need to include the derivative terms, as we have mentioned. 
Fine structures seen in the empirical density profiles, which come from the 
shell effects (see, e.g., a proton density dip at the center of a light ^^O nu- 
cleus), cannot be reproduced by the mean field approach. By imposing the 
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beta equilibrium, the most stable proton ratio can be obtained for a given 
mass number. Figure 1 (right panel) shows the mass- number dependence 
of the binding energy per nucleon and the proton ratio. We can see that the 
bulk properties of finite nuclei (density, binding energy and proton ratio) 
are satisfactorily reproduced for our present purpose. 

Note that we should adjust a slightly smaller value of the sigma mass 
than that one usually uses, i.e. 400 MeV, to get such a good fit. If we 
used the popular value of rn„ « 500 MeV, finite nuclei would be overbound 
by about 3 MeV per nucleon. Although the actual value of the; sigma 
mass (or the omega mass) has little relevance for infinite nuclear matter, 
it is important for finite nuclei and other non-uniform nucleon systems, 
since the meson mass characterizes the interaction range and consequently 
affects, e.g., the nuclear surface property. 




Figure 1. Left; the density profiles of typical nuclei. The proton densities (solid curves) 
axe compared with the experiment. Right: the binding energy per nucleon and the proton 
ratio of finite nuclei. 



4. Nuclear "Pasta" at Sub-nuclear Densities 

In the density region, where nuclei are about to dissolve into uniform nu- 
clear matter, it is expected that the energetically favorable mixed phase, 
which consists of nuclei and free nucleons, possesses interesting geometrical 
structures, such as rod-like and slab-like nuclei and rod-like and spherical 
bubbles, etc. These nuclei with exotic geometrical shapes are referred to as 
nuclear "pastas" . Phcnomcnologically the existence of the "pasta" phases 
instead of the crystalline lattice of nuclei would affect the mechanism of 
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supernova explosions and the glitch phenomena in pulsars. Due to these 
important consequences the "pasta" structure has been studied by number 
of authors. It is widely accepted that the "pasta" structure is realized 
basically due to the balance of the Coulomb energy and the surface ten- 
sion. However, the electron density has been always treated as an uniform 
background in the standard treatments. Here we study the nuclear "pasta" 
structure within our framework, which consistently treats the Coulomb po- 
tential and the electron distribution. 

4.1. Symmetric Nuclear Matter and the Coulomb Screening 
Effect 

First, we focus on symmetric nuclear matter (relevant to supernova matter 
at the initial stage of collapse) where the Coulomb screening effect by elec- 
trons is expected to be large. Figure 2 shows some typical density profiles 
in the Wigner-Seitz cells. The geometrical dimension of the cell is denoted 
as "3D" (three dimensional), etc. The horizontal axis in each panel denotes 
the radial distance from the center of the cell, and the boundary is indicated 
by the hatch. The nuclear "pasta" structures are clearly seen. Note that 
the electron density profile becomes no more uniform due to the Coulomb 
screening. 

The phase diagram of the matter structure is shown in Fig. 3 (left). 
The size of the cell i?ceii is optimized with precision of 1 fm, and the lowest 
energy configurations are chosen among various geometrical structures. In 
the figure, there never appears the spherical-hole configuration. This is one 
of the consequences of the Coulomb screening effect. It should be noted 
that appearance of the pasta structures is also sensitive to the choice of the 
effective interaction, as discussed by Oyamatsu et al.^^ 

To elucidate the Coulomb screening effect, there are two possible ways: 
one is to compare our results with those given by solving the equations of 
motion for fields and the density profiles neglecting the Coulomb potential 
Vcoui- Then the Coulomb energy is calculated by the use of the density 
profiles thus determined and finally added to the total energy, as in the 
simple bulk calculations. The optimum cell size is determined by this total 
energy including the Coulomb energy. The other way is to compare our 
results with those given by only discarding Vcoui in r.h.s. of the Poisson 
equation, while keeping it in other equations of motion. It is equivalent 
to the assumption of the uniform electron density distribution, which has 
been used in the previous studies;'^"® there protons interact with each other 
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Figure 2. Examples of the density profiles in the cell for symmetric nuclear matter 
(droplet, rod, slab, and tube). 
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Figure 3. Left: the binding energy per nucleon and the cell size of syiniiictric nuclear 
matter. Right: the same as the left panel with the uniform electron distribution. 
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and may form a non-uniform structure through the balance between the 
nuclear surface tension and the Coulomb interaction in a uniform electron 
background. Thus the density rearrangement effect is partially taken into 
account for protons, while it is completely neglected for electrons. The first 
way may be standard to extract the Coulomb effect and to compare our 
full calculation with the bulk one; actually this way has been taken in the 
context of kaon condensation^^ or hadron-quark deconfinement transition.^ 
However, we don't take the first way and dare to take the second way here 
to compare our results with the previous ones given by the uniform-electron 
calculation. 

We show in the right panel of Fig. 3 the results without the Coulomb 

screening (uniform-electron calculation). The region of each structure 
(droplet, rod, etc.) is different from that given by the full calculation. 
Especially, the "bubble" (spherical hole) appears in this case. Since the 
appearance of various geometrical structures and their region depend on 
the very subtle energy difference, the Coulomb screening has a significant 
influence on the sequence of the different pasta phases. 

Comparing the case of uniform matter with the case of the pasta phases, 
one can sec that the non-uniform structures reduce the energy. However, 
the Coulomb screening effect on the bulk EOS (difference between left and 
right panels of Fig. 3) is rather small. 

4.2. Nuclear Matter in Beta Equilibrium 

Next, we disciiss the nuclear matter in beta equilibrium, which is relevant 
to stable neutron stars. Figure 4 (left) shows the density profiles at several 
densities. Only the case for three dimensional (3D) geometrical structure 
is shown, since 2D and ID cases are energetically unfavored in our calcula- 
tion. One can sec the proton-enriched droplets embedded into the neutron 
sea at low densities. The EOS with the 3D phases is shown in the right 
panel of Fig. 4. In the beta equilibrium nuclear matter the effect of the 
non-uniform structure becomes much less compared to that for symmetric 
nuclear matter. Since the electron fraction is small, the Coulomb screening 
effect should not be remarkable. 

5. Summary and Concluding Remarks 

We have discussed the low-density nuclear matter structures "nuclear pas- 
tas" and elucidated the Coulomb screening effect. Using a self-consistent 
framework based on DFT and RMF, we took into account the Coulomb 



February 8, 2008 23:36 Proceedings Trim Size: 9in x 6in 



maxuyamaNM 



9 





0.08 




0.06 




0.04 


Q 


0.02 




0.00 
0.08 


[e 


0.06 




0.04 


Dens 


0.02 




0.00 
0.08 




0.06 


>i 


0.04 


Q 


0.02 




0.00 
0.08 


[b 


0.06 




0.04 


Dens 


0.02 
0.00 



~\ Pb=0.01 fin"' 3D 


1 


~^\^ PB=0.03fm-' 3D 
\ 

\ 


1 


\ 

\ 

\ 


^B=0.05 fm"' 3D 


Pb=0.08 fin-' 



^ 10 

> 



S 5 

I 



20 



ft: 10 



uniform assumption 

droplet 

uniform 



10 20 
Radius [fin] 



100 



0.05 



0.10 



Baryon density [fin ] 



Figure 4. Left: examples of the density profiles of nuclear matter in beta equilibrium. 
The electron density is too small to see clearly. Right: the binding energy per nucleon 
and the cell size. 



interaction in a proper way and numerically solved the coupled equations 
of motion to extract the density profiles of nucleons. 

First we have checked how realistic our framework is by calculating the 
bulk properties of finite nuclei as well as the saturation properties of nuclear 
matter, and found it can describe both features satisfactorily. 

In symmetric nuclear matter, we have observed the "nuclear pasta" 
structures with various geometries at sub-nuclcar densities. The appear- 
ance of the pasta structures significantly lowers the energy, i.e. softens the 
EOS, while the energy differences between various geometrical structures 
are rather small. So the Coulomb screening effect and the rearrangement 
of the charge density can affect such changes of geometrical structures in 
spite of that its absolute value is rather small in comparison with nuclear 
interaction energy. 

By comparing the results with and without the Coulomb screening, 
we have seen that the self-consistent inclusion of the Coulomb interaction 
changes the phase diagram. In particular one of the pasta configurations 
appears only when the Coulomb screening is switched off in our calculation. 
The effects of the Coulomb screening on the EOS, on the other hand, was 
found to be small. 
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We have also studied the structure of nuclear matter in beta equilibrium. 
There we have observed only proton-enriched droplet in the neutron sea. 
No other geometrical structures like rod, slab, etc. appeared. 

Detailed discussions about the Coulomb screening effect on the nuclear 
pasta phases will be reported elsewhere. 
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